Radiative corrections and parity nonconservation in heavy atoms 
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The self-energy and the vertex radiative corrections to the effect of parity nonconservation in 
heavy atoms are calculated analytically in orders Zo? and Z^o? \n{\c /ro), where Ac and ro being 
the Compton wavelength and the nuclear radius, respectively. The value of the radiative correction 
is —0.85% for Cs and —1.41% for Tl. Using these results we have performed analysis of the exper- 
imental data on atomic parity nonconservation. The obtained values of the nuclear weak charge, 
Qw = — 72.90(28)ei:p(35)theor for Cs, and Qw = —ll&.7{l.2)exp{iA)theor for Tl, agree with predic- 
tions of the standard model. As an application of our approach we have also calculated analytically 
dependence of the Lamb shift on the finite nuclear size. 
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Atomic parity nonconservation (PNC) has now been 
measured in bismuth ^ , lead , thallium ||^ , and cesium 
Analysis of the data provides an important test of 
the standard electroweak model and imposes constraints 
on new physics beyond the model, see Ref. [|. The 
analysis is based on the atomic many-body calculations 
for Tl, Pb, and Bi Q and for Cs 0, ||] (see also more 
recent Refs. 0)- Both the experimental and the 
theoretical accuracy is best for Cs. Therefore, this atom 
provides the most important information on the standard 
model in the low-energy sector. The analysis performed 
in Ref. Q has indicated a deviation of the measured 
weak charge value from that predicted by the standard 
model by 2.5 standard deviations a. 

In the many-body calculations 0, § the Coulomb in- 
teraction between electrons was taken into account, while 
the magnetic interaction was neglected. The contribution 
of the magnetic (Breit) electron-electron interaction was 
calculated in the papers [0 |l^. It proved to be much 
larger than a naive estimate, and it shifted the theoretical 
prediction for PNC in Cs. 

Radiative correction to the nuclear weak charge due 
to renormalization from the scale of the W-boson mass 
down to zero momentum has been calculated long time 
ago, see Refs. |l^, |lj]. This correction is always ac- 
counted in the analysis of data. However, another impor- 
tant class of radiative corrections was lost in the analysis 
of atomic PNC. This fact has been pointed out in the 
recent work [nsl that demonstrated that there are cor- 
rections ~ Zor caused by the strong electric field of the 
nucleus. Here Z is the nuclear charge and a is the fine 
structure constant. The simplest correction of this type 
is due to the Uehling potential. It has been calculated 
numerically in Ref. and analytically in our paper 



[ |l7[ . In the paper ||T^ we have also analyzed the gen- 
eral structure of the strong field radiative corrections. It 
has been shown that except of the usual perturbative pa- 
rameter Za there is an additional parameter ln(Ac/?'o) 
where Ac is the electron Compton wavelength and ro is 
the nuclear radius. 

In this letter we present the results of calculations of 
the strong field radiative corrections to the atomic PNC 
effect in orders Za^ and Z'^a^hi{Xc jr^). Details of cal- 
culations will be presented elsewhere |l^. Using our 
results we reanalize the experimental data. Agreement 
with the standard model is excellent. As an application 
of our approach we have also calculated the dependence 
of the Lamb shift on the finite nuclear size. Agreement 
of our analytical formula with results of previous compu- 
tations jl^, is perfect. 

The strong relativistic enhancement makes PNC ra- 
diative corrections different from previously considered 
radiative corrections to the hyperfine structure. The 
relativistic enhancement factor is proportional to i? ~ 
(Ac/(^aro)2(i-T), where 7 = - (zaf. The factor 
is R «3 for Cs and i? w 9 for Tl, Pb, and Bi The 
logarithmic enhancement of radiative corrections men- 
tioned above is closely related to the existence of the 
factor R. The Feynman diagram for the leading con- 
tribution to the PNC matrix element between pi/2 and 
Si/2 states as well as diagrams for radiative corrections 
are shown in Fig.^a) and Fig|^(b-f), respectively. The 
diagram Fig|^(b) corresponds to a modification of the 
electron wave function because of the vacuum polariza- 
tion. This correction calculated analytically in Ref. 
reads 
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FIG. 2: a{Za) self-energy and vertex radiative corrections. 
The zigzag line denotes Z-boson, and the dashed line denotes 
photon. 



FIG. 1; (a) Leading contribution to the PNC matrix element, 
(b-f) Radiative corrections. The double line is the exact elec- 
tron Green's function in the Coulomb field of the nucleus, the 
cross denotes the nucleus, the zigzag and the dashed lines de- 
note Z-boson and photon, respectively. 



where h = exp(l/(27) - C - 5/6), C « 0.577 is the Eu- 
ler constant, and / ~ 1 is some smooth function of Za 
independent of Tq. Hereafter we denote by 5 the relative 
value of the correction. So, Eq. (|^) represents the ratio 
of diagrams Fig.|l|(b) and Fig.|l|(a). 

The renormalization of the nuclear weak charge Qw 
from the scale of the W-boson mass down to g = was 
performed in Refs. [p^ [l4| . However, as it has been 
pointed out in Ref. ]17| |, atomic experiments correspond 
to g ~ 1/ro ~ 30MeV. The correction due to renormal- 
ization from g = to g = 1/ro is described by diagrams 
Fig.ll|(c) and Fig.ll|(d). It has the form O 



Scd = - 4sin2 Ow) ln(Ac/ro) ~ -0.1% , 



(2) 



where 9w is the Weinberg angle, sin^ 9w ~ 0.2230, see 
Ref. §. 

Diagrams Fig.|](c) and Fig.|l|(f) correspond to the con- 
tributions of the electron self-energy operator and the 
vertex operator, respectively. Each of these diagrams is 
not invariant with respect to the gauge transformation 
of the electromagnetic field while their sum is gauge in- 
variant. It has been demonstrated in Ref. ||l^ that the 
correction 5,, + 5f is of the form 



5ef^5e + 5f=A ln(&Ac/ro) + B 



(3) 



where A and B are some functions of Za^ and the con- 
stant b is defined after Eq.(|^) In the present paper we 
derive the functions A and B analytically in the leading 
approximation in the parameter Za, when A oc a{ZaY' 
and B cx a{Za). 

The simplest part of the work is calculation of the 
a{Za) contribution to B. It is convenient to use the 
Fried- Yennie gauge ||2^ together with the effective oper- 
ators approach |p3| where the corrections under discus- 
sion coincide with that for the forward scattering am- 
plitude, see Fig.||. The result of our calculation for the 
non-logarithmic term in Eq.(^) reads 



B = —a{Za) 
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For calculation of the function A (the logarithmic part 
in Eq. (||)) we have used the Feynman gauge. There 
are two contributions to A, the self-energy contribution, 
AsE, given by the diagram Fig.^e), and the vertex con- 
tribution Av given by Fig.|l|(f). It is convenient to repre- 
sent the self-energy operator as a series in powers of the 
Coulomb field of the nucleus, S Eq + Si -I- S2 + 
see Fig.^ and perform calculations in the momentum 
space. Distances -C r <C Ac where the logarithmic 
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FIG. 3: The electron self energy expanded m powers of the 
Coulomb field. The solid line is the free electron Green's func- 
tion, the cross denotes the nucleus, and the dashed line de- 
notes the photon. 

terms come from correspond to momenta m p <^ l/'"Oi 
where m is the electron mass. Therefore, the mass and 
the electron binding energy can be neglected in electron 
propagators. It is necessary to make the ultraviolet reg- 
ularization of the operators Sq and Si. Each of the op- 
erators depends on the parameter of the regularization. 
However, the contribution of Sqi = Sq + Si to Ase is 
independent of the regularization parameter because of 
the Ward identity. The operator S2 does not require any 
regularization. Our results for contributions of Sqi and 
S2 to Ase read 
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Next is the logarithmic contribution of the vertex op- 
erator described by Fig.|l|(f). The coordinate representa- 
tion is the most convenient for this part of the problem, 
and we use this representation. We have also used an 
integral representation of the Green function derived in 
Ref.Q (see also Ref. Q). Keeping in mind the contact 
nature of the PNC interaction and the logarithmic accu- 
racy of the calculation, one can demonstrate that only the 
angular momentum j ~ 1/2 is significant in the partial 
expansion of the electron Green's function in the exter- 
nal Coulomb field. A part of the vertex independent of 
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Za requires an ultraviolet regularization and hence is de- 
pendent of the regularization parameter. However, due 
to the Ward identity, the contribution corresponding to 
this part cancels out exactly the renormalization of the 
wave function, see e.g. Ref. |^. Another part is depen- 
dent on Za and gives the following contribution to the 
logarithmic term in Eq. (0) 
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Together with Eq. (||) this gives the final result for A 
a{Zaf /15 



A = -- 
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Thus, according to Eqs. (||), and (0) the total relative 
correction to the PNC matrix element due to the self- 
energy and the vertex operators reads 
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This correction is plotted in Fig. ^ versus the nuclear 
charge Z (the long-dashed line). The solid line in the 
same figure shows the total radiative correction to the 
PNC effect that includes both ^e/ (Eq. (||) and (5b (Eq. 

The radiative shift of the atomic energy levels (Lamb 
shift) depends on the finite nuclear size. This correction 
has a very similar structure to that of the PNC radia- 
tive correction since the effective sizes of the perturbation 
sources in both cases are much smaller than Ac . The self- 
energy and the vertex corrections to the finite-nuclear- 
size effect (SEVENS) for si/2-state have been calculated 
earlier analytically in order a{Za) in Refs.[p7|, |2^ . The 
corrections for lsi/2-, 2si/2-, and 2p]^/2"Sta'tes have been 
calculated numerically exactly in Za in Refs.|p^, 
The structure of higher in Za corrections and their loga- 
rithmic dependence on the nuclear size has not been un- 
derstood. We have applied our approach to the SEVENS 
problem and found the following expression for the Si/2~ 
state relative correction 
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Linear in Za term agrees with results of Refs. 
The logarithmic term coincides with that in Eq. (fel) for 
the PNC correction. Moreover, the logarithmic term in 
the SEVENS correction Ap for pi/2 state is also equal 
to that in Eqs.(|) and (|). The reason for this equality 
is very simple. The logarithmic terms come from small 



distances (r « Ac) where the electron mass can be 
neglected. When the mass is neglected the relative ma- 
trix elements for the PNC radiative correction and for 
SEVENS are equal. The correction A^ given by Eq. (||) 
is shown in EigJ by the dashed line. Results of the com- 
putations ||l^, |20| for Is and 2s states are shown by circles 
and diamonds, respectively. The agreement is perfect. It 
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FIG. 4: Relative radiative corrections (%) for the PNC and 
for the finite-nuclear- size effects versus the nuclear charge Z . 
The long-dashed line shows the correction S^f given by Eq. 
(Fig^(e,f)). The solid line shows the total radiative correc- 
tion to the PNC effect that includes both S^j and Sb (Eq.^, 
FigUfb)). The dashed line shows the correction As given by 
Eg. (m). Results of computations of As for Is and 2s states 
/ua, 1^3/ are shown by circles and diamonds, respectively. 

is interesting to find also the a(Za)-term in Ap. Using 
our result for the logarithmic term in Ap and results of 
numerical calculations of this correction for 2pi/2-state, 
Ref. |l^, we find that the coefficient in a{Za)-teTm 
is --1.9. 

It has been recently suggested in Ref. HI] that the 
following "precise relation" 5ef = (As -I- Ap)/2 is valid. 
Values of the corrections obtained in our work do not 
agree with this relation. The "derivation" in Ref. |2|] 
is based on the obviously wrong assumption that there 
is a gauge in which the vertex contributions to 6ef, As, 
and Ap vanish simultaneously. Although it is possible to 
vanish the vertex correction to each of these quantities by 
chousing an appropriate gauge, these are three different 
gauges. This is why the above relation is wrong. 

Now we can perform a consistent analysis of the ex- 
perimental data on atomic parity violation since all the 
contributions are known. In our analysis for Cs we have 
included the theoretical value of the PNC amplitude from 
Refs. !?[ H, ^, |l^ as well as —0.61% correction due to the 
Breit interaction iQ, —0.42% radiative correction calcu- 
lated in the present work, —0.2% neutron skin correction 
p2[ , —0.08% correction due to the renormalization of Qw 
from the atomic momentum transfer q 30MeV down to 
q = |17 , and -t-0.04% contribution from the electron- 
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electron weak interaction jlj. Using these theoretical 
results we obtain from the data |Q| the following value of 
the nuclear weak charge Qw at zero momentum transfer 



Cs: Qw ^ -72.90 ± (0.28)ex ± (0.35) 



theor • 



(10) 



This value agrees with prediction of the standard model, 
= -73.09 ± 0.03, see Ref. §. We have used the 
neutron skin correction in our analysis. However, in our 
opinion, status of the correction is not quite clear because 
data on the neutron distribution used in Ref. are not 
quite consistent with the data on neutron distributions 
obtained from proton scattering, see e. g. Ref. 

In the analysis for Tl we have included the theoreti- 
cal value of the PNC amplitude from Refs. as well 
as —0.88% correction due to the Breit interaction Q, 
—0.51% radiative correction calculated in the present 
work, —0.2% neutron skin correction, —0.08% correction 
due to the renormalization of Qw from the atomic mo- 
mentum transfer q ^ 30MeV down to q — Q ||l^, and 
-1-0.01% contribution from the electron-electron weak in- 
teraction |l^. Using these theoretical results we obtain 
from the data the following value of the nuclear weak 
charge at zero momentum transfer 



Tl: Qw' = -116.7 ±(1.2)e^± (3.4) 



theor • 



(11) 



This agrees with prediction of the standard model, Qw = 
-116.7 ±0.1, see Ref. §]. 

Concluding, we have calculated analytically for the 
first time all the strong electric field radiative corrections 
to the effect of atomic parity violation. This calculation 
has allowed us to perform a consistent analysis of the ex- 
perimental data. Agreement with the standard model is 
within 0.5(7. 

A.I.M gratefully acknowledge School of Physics at the 
University of New South Wales for warm hospitality and 
financial support during a visit. 



A.I.Milstein@inp.nsk.su 



[1] 



Email 
Email 
Email 

L. M. Barkov and M 



sushkov@phys.unsw.cdu.au 



I . S . Terekhov@inp.nsk.su 



S. Zolotorev, JETP Lett., 27, 357 
D. MacPherson et ai, Phys. Rev. Lett., 67, 



(1978); M. J 
2784 (1991). 

[2] D. M. Meekhof et ai, Phys. Rev. Lett., 71, 3442 (1993) 
[3] N. H. Edwards et ai, Phys. Rev. Lett., 74, 2654 (1995) 
P. A. Vetter et ai, Phys. Rev. Lett., 74, 2658 (1995). 



[4] C. S. Wood et ai, Science 275, 1759 (1997); S. C. Ben- 
nett and C. E. Wieman, Phys. Rev. Lett. 82, 2484 (1999). 

[5] D. E. Groom et al Euro. Phys. J. C 15, 1 (2000). 

[6] V. A. Dzuba, V. V. Flambaum, P. G. Silvestrov, and O. 
P. Sushkov, J. Phys. B 20, 3297 (1987); Europhys. Lett. 
7, 413 (1988). 

[7] V. A. Dzuba, V. V. Flambaum, and O. P. Sushkov, Phys. 

Lett. 141A, 147 (1989). 
[8] S. A. Blundell, J. Sapirstein, and W. R. Johnson, Phys. 

Rev. Lett. 65, 1411 (1990); Phys. Rev. D45, 1602 (1992). 
[9] M. G. Kozlov, S. G. Porsev, and I. L Tupitsyn, Phys. 

Rev. Lett. 86, 3260 (2001). 

V. A. Dzuba, V. V. Flambaum and J. S. M. Gingers, 



[10 

[11 

[12; 

[13 

[14 

[15 
[16 

[17 

[18 

[19 

[20 

[21 

[22; 

[23; 

[24; 

[25; 

[26' 
[27 
[28' 
[29 
[30 
[31 

[32' 
[33' 
[34 



hep-ph/0204131 

A. Dcrcvianko, Phys. Rev. Lett. 85, 1618 (2000) 
V. A. Dzuba, C. Harabati, W. R. Johnson, and M. S. 
Safronova, Phys. Rev. A 63, 044103 (2001). 
W. J. Marciano and A. Sirlin, Phys. Rev. D 27, 552 
(1983). 

W. J. Marciano and J. L. Rosner, Phys. Rev. Lett. 65, 
2963 (1990). 

O. P. Sushkov, Phys. Rev. A 63, 042504 (2001); 

W. R. Johnson, L Bednyakov and G. Soff, Phys. Rev. 

Lett. 87, 233001 (2001). 

A. L Milstein and O. P. Sushkov, Phys. Rev. A 66, 
022108 (2002). 

A. I. Milstein, O. P. Sushkov and L S. Terekhov, to be 
pubUshed. 

K. T. Cheng, W. R. Johnson, and J. Sapirstein, Phys. 
Rev. A 47, 1817 (1993). 

W. R. Johnson and G. Soff, At. Data Nuc. Data Tables 
33, 405 (1985). 

Parity nonconservation in atomic phenomena, 1. B. 

Khriplovich, Gordon and Breach, 1991. 

H. M. Fried and D. R. Yennie, Phys. Rev. 112, 1391 

(1958). 

M. I. Eides, H. Grotch, and V. A. Shelyuto, Phys. Rep. 
342, 63 (2001). 

A.L Milstein, V.M. Strakhovenko, Phys. Lett. A 90, 447 
(1982). 

R. N. Lee and A. L Milstein, Phys. Lett. A 189, 72 
(1994). 

K. Pachucki, Phys. Rev. A 54, 1994 (1996). 

K. Pachucki, Phys. Rev. A 48, 120 (1993). 

M. I. Eides , H. Grotch, Phys. Rev. A 56. R2507 (1997) . 



hep-ph/0206124 



M. Yu. Kuchiev and V. V. Flambaum 
S. A. Blundell, Phys. Rev. A 46, 3762 (1992) 
I. Lindgren, H. Persson, S. Salomonson, A. Ynnerman, 
Phys. Rev. A 47, 4555 (1993). 
A. Dcrcvianko, Phys. Rev. A 65, 012106 (2002); 
G. D. Alkhazov et al, Sov. J. Nucl. Phys. 26, 357 (1977). 
We found this value by rescaling the Breit correction for 
Cs |l^, see discussion in Refs. fisl , 0. 



